The Impact of Isospin Breaking on the Distribution of Transition Probabilities 
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In the present paper we investigate the effect of symmetry breaking in the statistical distributions 
of reduced transition amplitudes and reduced transition probabilities. These quantities are easier 
to access experimentally than the components of the eigenvectors and were measured by Adams et 
al. for the electromagnetic transitions in ^^Al. We focus on isospin symmetry breaking described by a 
matrix model where both, the Hamiltonian and the electromagnetic operator, break the symmetry. 
The results show that for partial isospin conservation, the statistical distribution of the reduced 
transition probability can considerably deviate from the Porter-Thomas distribution. 
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I. INTRODUCTION 



The spectral fluctuations in a rich variety of differ- 
ent physical systems show, if measured on the scale of 
the local mean level spacing, very similar features. This 
high degree of universality makes it possible to describe 
these fluctuations with random matrices. Random Ma- 
trix Theory (RMT) is a simple, schematic model in which 
the matrix elements of the Hamiltonian in some basis are 
replaced with random numbers. Apart from randomness, 
the only further input are the symmetries and invari- 
ances of the system, in particular time reversal invari- 
ance. It turns out that this assumption of full ergodicity 
or "chaoticity" leads, in many cases, to a complete and 
parameter free description of the spectral fluctuations. 
Such universal cases are said to be of Wigner-Dyson type. 
We refer the reader to the reviews in Refs. Orig- 
inally, Wigner had developed this approach in nuclear 
physics where it continues to find new applications. 

In recent years, interest has been focussed on devia- 
tions from the universal, parameter free Wigner-Dyson 
fluctuations. The deviations can have different reasons, 
such as regular effects competing with full chaoticity or 
the breaking of time reversal invariance, see Ref. 
Here, we wish to discuss deviations attributed to the 
breaking of isospin symmetry. We recall that symme- 
tries such as isospin or parity are, in contrast to time 
reversal invariance, associated with quantum numbers. 
In the context of spectral fluctuations and RMT, sym- 
metry breaking already received broad interest, see the 
compilation in Ref. g]. As early as 1960, Rosenzweig 
and Porter jsj analyzed atomic spectra by setting up a 
general random matrix model which describes crossover 
transitions between different angular momentum cou- 
pling schemes. In the late eighties, Mitchell et al. ^] 
measured and analyzed about 100 low lying states with 
known values of the isospin quantum number in the nu- 
cleus 26 Al. In Ref. §, this was discussed in the frame- 
work of a random matrix model which is a special case 
of the Rosenzweig-Porter model. An estimate for the 
statistical Coulomb matrix element, i.e. a measure for 



the degree of isospin breaking, could be obtained. Moti- 
vated by similar questions in molecular physics, Leitner 
and coworkers |^ performed a perturbative calculation 
of the spectral fluctuations in the random matrix model. 
More recently, new data were obtained in two statisti- 
cally highly significant experiments on the breaking of a 
point group symmetry in a resonating quartz block 
and the coupling of two chaotic microwave billiards 
Both cases, although physically very different, are statis- 
tically fully equivalent to symmetry breaking in quantum 
mechanics. Importantly, there is only one parameter en- 
tering the random matrix model. It is a unique measure 
e.g. for the root mean square statistical Coulomb matrix 
element. This illustrates that the random matrix model 
is the ideal tool to extract a root mean square symmetry 
breaking matrix element from the data. 

All these studies addressed the spectral fluctuation 
properties. Symmetry breaking, however, will also have 
an impact on the statistics of the wave Junctions and ob- 
servables sensitive to them. Recently, Adams, Mitchell 
and Shriner |^ collected reduced 7 ray transition prob- 
abilities from different experiments on ^^A\. As men- 
tioned above, this nucleus had already shown a strong 
deviation of the spectral fluctuations from the univer- 
sal Wigner-Dyson result due to isospin breaking [Q. 
The new results show that the distribution of the 
transition probabilities also considerably deviates from 
the Wigner-Dyson statistics, i.e. from the distribution 
that corresponds to full, parameter free "chaos". In the 
present contribution, we wish to discuss these results. 
To this end, we extend the random matrix model of 
Ref. 1^ to discuss transition probabilities. Similar in- 
vestigations were performed simultaneously and indepen- 
dently by Andersen, Ellegaard, Jackson and Schaadt px| ] 
for symmetry breaking in acoustic and elastomechanical 
systems. 

The article is organized as follows. In Sec. |l| we briefly 
review the experimental results on the reduced transition 
probabilities in ^^Al collected by Adams, Mitchell and 
Shriner [01. In Sec. [II, we discuss the random matrix 



model in the case that no symmetry is present, in par- 
ticular its predictions for transition probabilities. The 



1 



random matrix model for isospin breaking is numerically 



studied in Sec. IV. In Sec 



Dsp: 
ff 



the numerical simulation 



is approximated by a qualitative albeit analytical model. 
The data analysis is performed i n Se c. [v| . Summary and 
conclusion are presented in Sec. |VI]|. 



II. EXPERIMENTAL RESULTS 

Experimental reduced electromagnetic transition 
strengths between the excited states of the nucleus ^®A1 
have been collected by Adams et al. — see Ref. 
Their data involve levels between the ground state and 
the excitation energy of 8.067 MeV. In this region, states 
with isospin T = and T — 1 are found and isospin is 
known to be approximately conserved. 

The probability Bif of a transition from the initial con- 
figuration \i) to the final configuration |/) is the square 



B 



of the matrix element 



(1) 



(2) 



of the relevant transition operator O in a special basis. 

Approximately 180 levels and 1500 electromagnetic 
transitions are known. The fluctuations of the Bi f values 
shall be studied. To this end, their systematic depen- 
dence on the quantum numbers of the initial and final 
states must be removed. This has been done in Ref. ||] 
in the following way. The states in ^^Al have been char- 
acterized by their excitation energy E, spin J, parity 
TT, and isospin T. Transitions are characterized by their 
electromagnetic character X which may be 5 or , their 
multipolarity L and a label r which becomes isoscalar 
(IS) if AT = and isovectorial {IV) if AT = 1. Hence, 
both the transition operator O = 0{XL, t) and the tran- 
sition probabilities Bif = Bif[XL, r) are functions of X , 
L and t. We shall, however, not always write all these 
arguments. 

A transition sequence is defined as a set of reduced 
transition probabilities where the initial states have a 
common assignment J'^, T, the final states have a com- 
mon assignment J''^ , T', and the transitions have their 
three defining characteristics all the same. Thus the re- 
duced transition probabilities Bif — B{Ei, Ef) of a given 
transition sequence can be labeled by the energies Ei 
and Ef of the initial and final states, respectively. The 
aforementioned secular variations of the B values were 
removed by normalizing them to the local average value 
{B{Ei, Ef)) of B{Ei,Ef) so that the statistical variable 
used further on is 



viE^.Ef) 



B{E„E^ 



{B{E^,Ef))- 



(3) 



The local average is defined by help of weighting factors 
that are Gaussian functions of the excitation energies. In 



doing so one must, however, remove the systematic de- 
pendence of the local level distance D on the excitation 
energy; i.e. one measures the energy in units of D and 
works with a dimensionless energy 



(4) 



The details are given in Ref. [|ri|. The local average of 
the B values of a given transition sequence is then 



{B{e^,ef)) = 



Ee.s'B{e,e') e 



V8 



(5) 



The fact that the variance of the Gaussian functions has 
been chosen equal to 4 is discussed in Ref. [pi]] . 

It is clear that the definition of the quantities y = 
y{Ei,Ef) requires the necessary spectroscopic informa- 
tion for the relevant states in ^^Al. Furthermore, the 
local averages D and {B{Ei, Ef)) require that there is 
a minimum number of members in the sequences from 
which these averages are derived. These restrictions have 
finally led in Ref. to an ensemble of 873 values of y. 

The authors of Ref. have found it convenient to 

transform y to the logarithmic variable 



logio y 



(6) 



The experimental distribution of z is given by the his- 
togram on Fig. Although the widths Az/j, of the bins 
fc = 1 • • • 18 in Fig. |^ vary, the ordinate is the probability 
density with respect to z: The height pk of the k-th bin 



Pk 



Azk 



(7) 



where Nk is the number of cases falling into the fc— th 
bin. 

Since the states \i) and |/) are believed to be "very 
complicated" , it is natural to assume that the matrix ele- 
ments (0) have a Gaussian distribution; see the discussion 



in Sec. |IIB. This entails that y has a Porter-Thomas 



distribution with unit average value, i.e. 



1 exp(-y/2) 



Transforming this to the variable z yields 
In 10 



Piz) 



10^/2 exp(-1072). 



(8) 



(9) 



A comment on the notation is in order: For every prob- 
ability density we shall use the symbol P. The fact that 
we deal with different functions will be clearly indicated 
by the argument of the probability density in question. 

The distribution (|^) is given on Fig. ^ for compari- 
son with the data. The distribution (pl) is normalized 



2 



to unity; the histogram of the data, however, has been 
normahzed to 



18 



^PfeAzfc = 0.83 



(10) 



k=l 



by the authors of Ref . Q . This is motivated ||T^ by the 
fact that / dz P{z) over the range of the data amounts 
to 0.83 and approximately takes care of the upper and 
lower detection thresholds. 

The point of Ref. is the shift of the experimental 
distribution with respect to the Porter-Thomas distribu- 
tion P{z). The maximum of the histogram occurs around 
z = —0.5. The maximum of P{z) occurs at z = 0. The 
authors of Ref. |^ have conjectured that the discrepancy 
between the experimental data and the Porter-Thomas 
distribution may be a consequence of isospin breaking in 
^^Al. The following sections of the present paper focus 
on this conjecture. In the next section, a random ma- 
trix model of isospin violation — and more generally of 
symmetry breaking — is presented. 
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FIG. 1. The distribution of experimental reduced transi- 
tion probabilities B{XL,t) in ^"^Al from Ref.[9]. The quan- 
tity B has been transformed to the logarithmic variable z of 
Eq. (6). The curve labeled Porter- Thomas is the distribu- 
tion (9). 



III. 



RANDOM MATRIX MODEL - 
OF NO SYMMETRY 



THE CASE 



In Sec. [II A, we briefly review the distributions of 
wave functions and widths for the standard Wigner- 
Dyson case in which no symmetry is present. A more de- 
tailed presentation with further references can be found 
in Ref. H. In Sec. [II B, we work out the distribution of 



A. Wave functions and decay amplitudes 

If the system is invariant under time reversal, the wave 
functions can be chosen real and the N x N random ma- 
trix H modeling the Hamiltonian is real and symmetric. 
The matrix elements are Gaussian distributed random 
numbers and H is said to be in the Gaussian Orthogo- 
nal Ensemble (GOE). In the physically relevant limit of 
infinitely many states, N ~f oo, the fluctuation proper- 
ties are of Wigner-Dyson type. A wave function is mod- 
eled by an eigenvector Ui, i = 1- ■ ■ N oi H, i.e. we have 
Hui — EiUi where Ei is the eigenenergy. 

We are interested in the probability density P{a) of 
finding the value a with — 1 < a < +1 for a component 
Uim of the eigenvector Ui. For finite N one finds 



r(A^/2) 



0Fr((A^-l)/2) 



(l-a2)(^-3)/2 



The second moment of this distribution is 



(11) 



For a large number of levels one obtains a Gaussian with 
variance 1/N, i.e. 



■ exp 



27r 

G{a,N-^'^) . 



N . 



(12) 



Here, we have introduced the notation G(a, a) for a Gaus- 
sian with variance a^. 

The wave functions are rarely accessible in an exper- 
iment. However, other observables, such as scattering 
matrix elements or partial widths, sensitively depend on 
them. Consider the scattering from a state i with wave 
function Ui into a channel c with channel wave function 
Xc- The corresponding reduced partial width amplitude 
7ic can be written as 



N 
m — 1 



(13) 



where J^c is the overlap integral between the m— th 
canonical basis vector and Xc- Again, the probability 
density can be worked out 



T{N/2) 



[TTN-i^^^'"^ T{{N -l)/2) 



1 = 



{N-3)/2 



(14) 



the transition matrix elements for the same case. 



where we suppress the indices i and c. The second mo- 
ment reads 7^ = 7^, = N^^ J2m=i ^mc- We notice that 
the functional form of this distribution agrees with that 
of Eq. (pd]). Thus, in the limit of large N, one again finds 
a Gaussian. Usually, one introduces the partial width 
r = 7^ which can be measured. The relative partial 
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width y = T /T with F = 7^ is distributed according to 
the Porter-Thomas law (||). 

We notice that the Porter-Thomas law or, equiva- 
lently, the Gaussian for the partial width amplitudes 



results from a large N expansion of the distribution (14). 
Alternatively, one may derive these large N results by 
using the central limit theorem: The partial width ampli- 
tudes are, according to Eq. (|l^), given as a linear combi- 
nation of N components ■ If — as assumed in — 
all these components are independently distributed, the 
distribution of the partial width amplitudes approaches 
a Gaussian for large N. This line of arguing does not 
use the fact that the distribution of every single compo- 
nent Uim is Gaussian for large N. It would apply for any 
smooth distribution of Uim, provided it does not sensi- 
tively depend on A^. This subtle point will be important 
in Sec. ffVAl. 



B. Transition matrix elements 



The results compiled in Sec. |III A| apply to the par- 
tial widths and to their amplitudes. However, in the ex- 
periments on ^^Al, electromagnetic transition probabili- 
ties Bif(XL) were measured which are squares of tran- 
sition matrix elements according to Eqs. (l|) and In 
Rcf. it is argued that the distribution of the tran- 

sition matrix elements Wif is, once more, of the form 
and (p^. In the present section, however, we give a 
derivation, valid for large N, which is well suited for the 
discussion of isospin breaking in Sec. [V. 

In our model, the initial state \i) and the final state |/), 
are represented by the eigenvectors Ui and Uf, respec- 
tively. The electromagnetic transition operator 0{XL) 
is modeled by one fixed random matrix. It is not nec- 
essary to consider an ensemble of such operators. Thus, 
the transition matrix elements Wif (XL) are given by 



W,f{XL) = uJO{XL)ui 



\ufnOnm{.XL)u 

V 



(15) 



The transition probabilities read Bif{XL) — 
\Wif{XL)\'^. Although these quantities differ from the 
partial widths and their amplitudes, it is intuitively ob- 
vious that, for large A'', the Wif {XL) are Gaussian dis- 
tributed and that the Bif{XL) obey a Porter-Thomas 
law: Since we always consider i ^ f, the Wif (XL) are 
linear combinations of (products of two) independent 
variables and, therefore, the central limit theorem ap- 
plies. 

More precisely, the distribution of the matrix elements 
Wif reads 



P{W.f) = / d[u,]P{u,) J d[uf]P{uf) 
5{W,f^uJOu,) 



where we suppress the argument XL of Wi / and O. Since 
we may assume that A^ is large, the distributions of the 
components Uim and Ufm take the form 



P(u) 




exp 



m— 1 



N 



■N 



exp 



A^ 



(17) 



with u standing for either Ui or Uf. We notice that the in- 
tegration domain for each of the eigenvector components 
is the interval [—1,-1-1]. However, since we consider the 
large A" limit, the distributions (^ are so sharply peaked 
at the origin that we may extend the domain of integra- 
tion to the entire real axis. Using the Fourier transform 
of the S function in Eq. (|l^) we have 



PiW^f) = ^ r dte^p{itW,f)Rit,0) . 



(18) 



The Gaussian integrals absorbed in R{t, O) can be done 
easily. 



I 2N 



exp 



-yK Uf] 



In itO/N 
itO/N In 

det'^^^ {In + t'^0^/N^) 



Uf 



exp 



(^~Ur\Yi{tN + t^0^lN'')^ . (19) 



We are interested in large A^. This allows us to expand 
the logarithm and to keep only the first term which is 
t'^O^/N^. Thus we have for large A^, 



R{t, O) = exp 



Q2 



(20) 



where the second moment of the operator is defined by 

(21) 



0- = ^trO^ 



We assume that the matrix elements Onm{XL) do not 
depend on A^. Hence, collecting everything, the Fourier 
transform in Eq. ( |l^ ) yields 



P{W,f)^G{W,f,0^ ) , 
which is the expected result. 



(22) 



IV. RANDOM MATRIX MODEL — THE CASE 
OF ISOSPIN BREAKING 



In Sec. IVA, the model is set up for the case that 



(16) isospin is partially conserved. In Sec. IV B, numerical 
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simulations are described. Note that we speak of isospin 
breaking because we have in mind the experiment of 
Sec. ||. The model, however, applies to any other quan- 
tum number that is partially conserved. 



A. Definition of the model 

We consider two isospin values T = and T — 1. If 
isospin were fully conserved, the Hamiltonian H would 
be block diagonal. The Coulomb interaction, however, 
destroys this symmetry. In Ref . ||^ , this was modeled by 
using random matrices of the form 



H = 



H{0) 













He' 


+ a 


HE 






(23) 



where H(j), j = 0, 1 are independent GOE matrices with 
dimension Nj . The total level number is iV = iVo + iVi . 
The Nq X iVi matrix He accounts for the Coulomb in- 
teraction. It is real without further symmetries and has 
Gaussian distributed entries. In this model, the param- 
eter a is proportional to the root mean square Coulomb 
matrix element. Since H{j) and He have positive and 
negative entries with equal weight, all observables can 
only depend on the modulus of a, but not on its sign. 
Thus, in the sequel, we restrict ourselves to a > 0. We 
recall that the spectral fluctuations are measured on the 
scale of the local mean level spacing D. Thus, the rele- 
vant parameter governing the spectral fluctuations is 



A = a/D. 



(24) 



For a = we have a non-interacting superposition of 
two independent GOE's. In this case, the distribution of 
the eigenvector components of the full matrix H, i.e. in- 
cluding the many exact zeros, is given by 

PnqNi (a, 0) = glPNo (a) + glPm («) + 25051-^(0) , (25) 

where we have introduced the fractional level numbers 

go = No/N and gi ^ N^/N . (26) 

The distributions P^v (a), j = 0, 1 have to be taken as 
either Eq. ( pi] ) or (p^). Obviously, the total distribution 
PnqNi (a, 0) is properly normalized. 

As discussed in Ref. Ql, the variances of the distribu- 
tions for the matrix elements are chosen in such a way 
that the spectra of H{Q) and H{1) have the same length 
and that H becomes a full N x N GOE matrix [|l5| for 
a = 1 if A'o = A^i . This also means that 



PN/2,N/2{a, 1) = P7v(a) 



(27) 



For a = 1 but Nq ^ N\ there are some deviations from 
the pure GOE results. 

For arbitrary a, the distribution PNgN-^{a, a) is not 
known analytically. A qualitative model is presented in 
Sec. (see Eq. (|^)). 



The transition matrix elements Wif{XL) and the cor- 
responding transition probabilities Bif{XL) depend not 
only on the eigenstates Ui of H, but also on the transition 
operator 0{XL). Since the latter contains the effective 
charges or the magnetic g factors of proton and neutron, 
it causes isospin breaking in addition to the isospin break- 
ing built into ( p3). Thus, in the same isospin basis that 
was used in Eq. (23|), we may model the transition oper- 
ator by a matrix of the form 



0(0) 

0(1) 



Oe 
01 



(28) 



Here, each of the matrices 0{j), Oe will be modeled by 
one fixed random matrix — as described in Sec. [VB. 
For later convenience, O has been written as a function 
of the two parameters /?/ and (ie- However, only the 
ratio 



P = Pel Pi 



(29) 



is important since the total strength of O drops out of 
the observable y of Eq. (^). For /3 = 0, configurations 
with different isospin values are not mixed by O while 
for P ~ 1 the mixing is maximal. For /? — s- 00 the opera- 
tor O couples configurations with different isospin only. 

We are interested in the distribution of the Wif{XL) 
and the Bi AXL). According to the discussion in the pre- 
vious Sec. In, one is tempted to argue as follows: Since 
Eq. (^5|) is general and also holds in the presence of con- 
served or broken isospin, the transition matrix elements 
Wif{XL) are always a linear combination of products 
of independently distributed sets of variables and 
u fn- Thus, the central limit theorem should apply and we 
should obtain a Gaussian distribution for the Wif{XL). 
This should even be true for any smooth distribution for 
the wave function components and u/„. In other 
words, the specific form of these latter distributions, 
which are non-trivial functions of the mixing parame- 
ter a, is not important and we would always expect a 
Gaussian distribution for the transition matrix elements 
Wif{XL). However, this reasoning is incorrect, because 
it does not make use of a specific basis. In other words: 
The result of this reasoning is an orthogonally invariant 
distribution. Within the random matrix model ( p3| ) one 
sees that this cannot be true. The statistical properties 
of the Hamiltonian H are not orthogonally invariant — 
except for the special case of no isospin symmetry at all. 
One sees the lack of orthogonal invariance very clearly in 
the distribution (|2^) of the eigenvector components in the 
case a = 0: The (5-distribution will appear only if the ba- 
sis vectors have well defined isospin. In general, the wave 
function components are functions of a. These functions 
parametrically depend on N. This A^-dependence com- 
petes with a large N expansion needed in the derivation 
of the central limit theorem — whose premises are there- 
fore violated. We illustrate this in Sec. IV B by numerical 
examples. 
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B. Numerical simulation 

The distributions of the reduced transition amphtudes 
and reduced transition probabihties were numericaUy in- 
vestigated. To this end, random Hamiltonians with the 
structure of Eq. ( p3| ) have been constructed. The dimen- 
sions A'o and A^i were chosen to be 100 so that H has 
the dimension N = 200. The elements of H{j), j = 0, 1, 
were selected by a generator of Gaussian random num- 
bers such that the second moments were 



{Hc)fj.miHc) 



(30) 



This ensures that _ff is a full GOE matrix for a = 1. The 
mean spacing Dj of the eigenvalues of H{j) is then (in 
the center of the spectrum) 



-1/2 



= 0.314. 

The mean level spacing D oi H is 

D = {D^^ +D^^)-'^ = 0.157. 



(31) 
(32) 



This value changes little, when a is varied between and 
0.157. This range of values was considered in the present 
numerical simulations. The parameter A of Eq. ( p^ ) then 
covers the range of < A < 1.0. Ten Hamiltonians were 
constructed and diagonalized for each value of A. 

The transition operator O was constructed very much 
as the Hamiltonians; i.e. its elements were selected by 
the random number generator. The variances were cho- 
sen in complete analogy with Eqs. (pO|). However, for 
each value of /3, one operator has been generated. The 
range of < /3 < 1.0 was covered. 

For given parameters A, (3 the transition ampli- 
tudes ( |l5| ) have been calculated with the indices i, f run- 
ning over the eigenvectors of the ten Hamiltonians. The 
set of 2 • 10^ numbers Wif with i > f forms the numeri- 
cally generated data. 

In a first step, we have checked whether these data 
follow a Gaussian distribution. To this end one can bin 
the set of Wif and compare the histogram to a Gaussian; 
similarly one can compare the distribution of the variable 
y of Eq. (m to the distribution (|q) or the distribution of 
z of Eq. (o) to P{z) given in (PJ! In order to make the 
comparison — as much as possible — independent of the 
parametrization of the statistical variable, the general- 
ized entropy 



S 



dx p{x) In 



P{x) 



(33) 



has been used to express the difference between the distri- 
bution p of the data and the ansatz P. This expression is 
independent of the parametrization x. Since both, p and 
P, are normalized, S is never positive. It vanishes if and 
only iip= P. The data do of course not provide a con- 
tinuous probability density p(x) but rather a histogram 



where the fc— th bin is centered at Xk and contributes 
the probability p(xk)l^Xk. Therefore they provide the 
approximation 
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^ AxkP(Xk) In— r — r- 

P{Xk) 



2^ 

k 



(34) 



to expression (p3D. To the extent that (^4|) is a good 
approximation to (^) , expression (^J) is independent of 
the parametrization. This holds in the present case to 
within a few times 10""* due to the large ensemble of 
numerical data. In Tab. ^, the entropy (Q) is given for 
a few values of (3 and A in the range of < /9, A < 1. 
One finds that the distribution of the data does depend 
on the parameters A and /3 which govern the symmetry 
breaking. The entropy approaches zero, when /? or A ap- 
proach unity, i.e. when the symmetry is strongly broken. 
Hence, in this case the transition amplitudes Wif have 
a Gaussian distribution. This does not hold if (3 and A 
have values much smaller than 1, in which case neither 
the Hamiltonian ( p3| ) nor the operator (|8|) significantly 
breaks the symmetry. 

TABLE I. Entropy for diflerent values of the parameters 13 
and A. The parameter a is also shown. 



A 


13^ 
ai 





0.02 


0.06 


0.1 


0.4 


1.0 








-1.05 


-0.75 


-0.53 


-0.37 


-0.04 


-0.0003 


0.012 


0.002 


-0.54 


-0.53 


-0.40 


-0.30 


-0.04 


-0.0003 


0.019 


0.003 


-0.43 


-0.42 


-0.34 


-0.26 


-0.04 


-0.0003 


0.031 


0.005 


-0.35 


-0.34 


-0.29 


-0.23 


-0.03 


-0.0003 


0.125 


0.02 


-0.06 


-0.06 


-0.06 


-0.06 


-0.01 


-0.0003 


1.0 


0.16 


-0.0006 


-0.0006 


-0.0006 


-0.0004 


-0.0004 


-0.0004 



The figures |2| and |3| illustrate in which way the dis- 
tribution of Wif deviates from a Gaussian. Consider the 
case oi (3 = = A, i.e. the absence of symmetry breaking, 
displayed in Fig. ||A. There is no transition between dif- 
ferent symmetry classes and therefore about one half of 
the matrix elements Wif vanishes. This leads to the nar- 
row peak in Fig. |^A comprising one bin centered at zero. 
This peak is superimposed over the distribution of the 
matrix elements connecting states in one and the same 
symmetry class. Fig. |^A thus corresponds to Eq. ( [25| ) 
with Nq — Ni. When the symmetry is completely bro- 
ken — in Fig. |D — the matrix elements Wif have a 
Gaussian distribution. Intermediate situations are dis- 
played in the parts B and C of Fig. ^ 

The intermediate situation is especially well visualized 
on Fig. |[ where the distributions of Fig. || have been 
transformed to the logarithmic variable 



log 



10 



W^f 



(35) 



which was introduced in Eq. 
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FIG. 2. Numerically generated distributions of transition matrix elements Wif for four different pairs of mixing parameters 
(3 and A. The matrix elements are normalized to their root mean square value. The binning of the data is the same in the four 
parts of the figure. The scale of the ordinate is different in the upper and lower parts. 



Here, the distribution displays two peaks. One of them 
— centered at z = — stems from the matrix elements 
connecting states within one symmetry class. The other 
one — centered at negative values of z — is due to the 
matrix elements connecting states from different symme- 
try classes. In case of complete symmetry breaking the 
two peaks merge as in Fig. ^D. This distribution corre- 
sponds to Eq. (^. In the present parametrization the 
case of no symmetry breaking (Fig. ||A) cannot be dis- 
tinguished from the case of complete symmetry breaking 
(Fig. |d). This occurs because z —oo when W ^ 0: 
hence, the matrix elements with the value W — are 
found nowhere on the z— axis. 

In Sec. |l^ it has been described how the secular vari- 
ations of the statistical properties of the experimental 
data of Ref. ||] have been removed by help of the local 
average in Eq. (^) . We have applied the same procedure 
to some of the present numerically generated data sets. 
Introducing a local average {Bif) instead of the global 



average W^. 



'if — Bif used in Figs. does not change 
the present results; i.e. along the spectrum of eigenstates 
there are no secular variations of the statistical properties 
of the eigenfunctions. 

The numerical simulations clearly demonstrate that 
there indeed is an effect of symmetry breaking on the 
distributions. Hence, it is shown that the complexity of 
the wave functions alone is not enough to ensure Gaus- 
sian distributions of wave function, width and transition 
amplitudes, as a naive application of the central limit 
theorem would imply. 



V. QUALITATIVE EVALUATION OF THE 
RANDOM MATRIX MODEL WITH BROKEN 
SYMMETRY 



In the present section, we discuss the random matrix 
model in a heuristic fashion. This is an extension of the 
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Sec. II] to the case of broken symmetry. Although we 
beHeve that the model of Eqs. (|2^ , ^8| ) can be completely 
solved analytically, we refrain from tackling that task. As 
we shall see, qualitative reasoning will provide sufficient 
understanding of the statistical behavior of the obscrv- 
ables. Thus, our analytical derivations will be approxi- 
mate and only exact in some limiting cases. Of course, 
we shall compare our approximate results with the above 
numerical simulations of the full model. 

To obtain the distribution P{W) we use the ansatz ( p^ ) 
for the transition operator and we need an ansatz for the 
distributions of the wave functions. The latter ones — in 
our model of the eigenvector components — are the en- 
tries of the orthogonal matrix U which diagonalizes the 
Hamiltonian, 



U = 



U{0) C/(01) 

um u{i) 



(36) 



where the matrices U{j), j — 0,1 are Nj x Nj orthog- 
onal and U{01) and U{10) are Nq x A^i and A^i x Nq 
matrices, respectively, with the condition that the to- 



tal U is orthogonal. Motivated by the numerical sim- 
ulations in the previous section and by the general ex- 
perience with wave function distributions, we assume 
that the distributions of the entries can be approximated 
by Gaussians. More precisely, we estimate that the 
distributions of the U (j) matrix elements are roughly 
given by G{a, kjN''-^/^), j = 0,1, while those of the 
U{01) and U{10) matrix elements are roughly given by 
G{a,r]N^^^'^). Thus, we reduce the problem to the de- 
termination of the three parameters Kq, ki and rj. It is 
reasonable to assume that the former two are functions of 
the latter one, i.e. we set j{ri). This is motivated 

by the fact that ij can replace the mixing parameter a 
(or A) and that — apart from the level numbers Nq, Ni 
— the distribution of the elements of [/ is a function of 
the symmetry breaking. The final task is then to find the 
relation between the purely phenomenological parameter 
ry and the parameter a which measures the root mean 
square symmetry breaking matrix element and thus has 
a direct physical interpretation. In other words, we even- 
tually have to specify the function rj = r]{a). 
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FIG. 3. Distributions of numerical transition probabilities Bif. These are the same data as in Fig. 2, but reparametrized in 
terms of the variable z of Eq. (6). 
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To construct a consistent guess for the functions Kj = 
Kj{r]), we notice that, in line with our basic assumptions, 
the combined distribution for the matrix elements of the 
total U reads 



PNoNiia,ri) = g^G a. 



N 

2.9ogiG ( a. 



giG [ a. 
V 



(37) 



This form is consistent with the exact results (g^) for ar- 
bitrary No and iVi and (jH) for Nq — Ni. We mention in 
passing that we always deal with Gaussian approxima- 
tions of distributions which, in some limiting cases, are 
for finite level number given by Eq. (pT|). 

Let us consider the case a = which yields the distri- 
bution (Ea). Comparison with Eq. (B7h shows that 



77^(0) = 



(38) 



and 



4(0) 



In the case of a 
leads to 



and 



- l/gj , J =0,1. (39) 
1 and iVo = A^i , the exact result ^) 



ry2(l) 



1 



(40) 



(41) 



Kj(l) = l, j = 0,l. 

Moreover, the functions Kj = Hj{ri) ought to be even in 
ry. This suggests the ansatz 



1 

9j 



1 



1 

9j 



V 



J = 0, 1 . 



(42) 



It is reassuring that similar combinations of parameters 
do also show up in exact calculations of the spectral cor- 
relators. 

A reasonable form for the function rj = ri{a) remains 
to be given. The function will parametrically depend on 
go and gi . We construct it by comparing the ansatz ( |37| ) 
with the distribution of the coefficients of the eigen- 



IV B 



functions Ui from the numerical simulation of Sec 
For each value of A = a/D, the optimum 77 was deter- 
mined by maximizing the entropy (^). We have not at- 
tempted to assign errors to the values of 77 thus found. In- 
stead we have convinced ourselves by inspection that the 
ansatz ( p7| ) reproduces the numerical distribution reason- 
ably well. A typical example is given by Fig. ^. In this 
way the function 77 — ??(q^) was generated numerically. 
For the case of A'o = Ni is is displayed on Fig. || together 
with the analytical expression 



772(a) = l-exp(-a/0.157) 



(43) 



that well approximates it. The numerical constant in the 
exponential is here equal to the level distance of Eq. (32). 



We emphasize, however, that the relation ( p3| ) is purely 
phenomenological. In Sec. we shall see that (^ ) can 
successfully be applied to the case of ^^Al, where D is 
different from (|32|), see Eq. (^) below, but the differ- 
ence is not too large. Thus it seems that the numerical 
constant is indeed related to the level distance. This 
quantity shows a significant N dependence, see Eq. (|3l]). 
It is at this point where the premises of the central limit 
theorem are violated for the calculations of the distribu- 
tions of the transition matrix elements. They cannot be 
a single Gaussian for all values of the parameters. 
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No=Ni=100 
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FIG. 4. The distribution of the components o of the eigen- 
functions of H. The parameters of the numerical simulation 
are given on the figure. This yields the histogram. The curve 
is the function (37) with the optimum -q. 

Having fixed all parameters in our qualitative model 
for the distribution of the eigenvector components, we 
can now calculate the distributions for the transition ma- 
trix elements Wif = uJOui. We emphasize that this 
involves, apart from the large TV limit, no further ap- 
proximation. We consider i ^ f, such that the initial 
and final states Ui and Uf are always different columns 
of the matrix U in Eq. (pq). Due to isospin breaking, we 
have to distinguish four structurally different possibili- 
ties: in case (MOO) both, Ui and Uf, are among the first 
A'o columns. We write 



and Uf 



MO) 
«/(10) 



(MOO) 



In case (Mil), both are from the last iVi columns, 



U^{^) 



and 



^■(01) 



(Mil) 



In the cases (MOl) and (MIO), Ui and Uf are from differ- 
ent parts in U, we have 
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and 



U^{10) 



U^{0^) 
U^il) 



and Uf 



and 



u/(01) 



Uf {10) 



(MOl) 



(MIO) 



Thus, for the case {Mjj'), j, j' = 0, 1 we will obtain a dis- 
tribution Pjji{Wif). Because of the symmetries O-^ — O 
and Wfi =Wif, the cases (MOl) and (MIO) must yield 
identical distributions. 
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FIG. 5. The relation between a and rj for TVo = TVi. The 
crosses are from the numerical simulation. The curve is a 
phenomenological representation of the numerical results. 

We proceed b y gene ralizing the evaluation of the inte- 
gral ( |l6| ) in Sec. Ill B| . Again, the distribution 



Pn'^f)^ J d[u,]P{u,) J d[uf]Piuf) 
is rewritten in the form 



Pif (W^f) = ^ r dt exp (ttW.f) {t, O), 



(44) 
(45) 



where the main difficulty is the calculation of the function 

Rjrit,0)= f d[u,]Piu,) f d[uf]Piuf) 



exp {-ituJOu,) . (46) 

In the case (MOO), the distributions of the eigenvectors 
are according to our ansatz 



I N 



-No 



( -7r-2"r(o)u»(o) 



/ TV 
2^ 



-Ni 



exp 



-^wf (10)u.(10) ) (47) 



for the initial state and 

I No 



P{uf) 



I N 
2'kkq 



exp 



-Ni 



N 



«/(0)"/(0) 



N 



N 



exp -— u}(10)u/(10) (48) 



for the final state. The case (Mil) is treated similarly. 
The case (MOl) has a different structure: The initial state 
distribution is as in Eq. ( ^ ) , but the final state distribu- 
tion reads 



-No 



N 



e^pl^ uj{l)uf{l) . (49) 



N 



2'kkI \ 2k{ 



However, one sees that in every case we may cast Eq. 
into the form 



i?,y(t,0) = Qy / d[u,] j d[uf] 

N 



exp 



Ui 



= C]y dct-1/%,,. 



(50) 



Here, Cj^' and C'^y collect all normalization constants. 

It is straightforward to write down the matrices ¥ijy . 
We define the diagonal matrices 



Do = diag ( ^Iato, ^Iwi 

Di = diag ( ^Iato, ^Ia^i 
.^7 '^i 



and find 



Dj itO/N 
itO/N Dy 



(51) 



(52) 



in the four cases. An easy calculation yields 
Rjf{t,0) = det"^/^ In + —Di^ODr^O 



iV2 J 



1 



exp ( --tr In I Iat + —Dj'OD-'O ] ] (53) 



which leads in the large N limit to 



Rj-j' {t, O) = exp 



CI. 



n ^2 



(54) 
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The effective second moment O^^, of the operator O is 
defined by 



(55) 



Hence, by performing the Fourier transform in Eq. (45) 
we arrive at 



P33'iWrf) = G{W,f,Ol,^'^) , 



30 



i.e. at a Gaussian in all four cases. The Gaussian form 
does not come as a surprise. The crucial and non-trivial 
result of this calculation is the explicit form of the effec- 
tive second moments (|55|). Denoting the second moments 
of the blocks in the operator (^ by 



1 



1 

N2- 



(57) 



we find after a simple computation in the four cases 



(58) 



The distributions Pjf{Wif) of Eq. (|5q) describe the 
cases {Mjj'). These were defined by assuming that the 
structure of the initial and final states, i.e. their place of 
origin in the matrix ( p6| ) is known. In a set of experi- 
mental data, this information is usually lacking. Thus, 
the full distribution must be a proper superposition of 
the functions Pjj'{Wif). The weights in this superpo- 
sition can be easily constructed: In the case (MOO), 
there are A'o states and therefore Nq{Nq — 1) transi- 
tions. We recall that we exclude recoupling transitions 
with i — f. Since, altogether, there are N states and 
N{N — 1) transitions, the weight of the case (MOO) is 
NoiNo - l)/^(^ - 1), which yields for large level 
numbers. In the case (Mil), one obviously finds gf. Sim- 
ilarly, we find NoNi/N{N — 1), i.e. gogi for large level 
numbers, in both of the cases (MOl) and (MIO). Thus, 
collecting everything we eventually end up with 



PiW.f)=glG{W.f,0, 



-25051 G(M^./, 02/^'), 



(59) 



where the individual variances are explicitly given in 
Eq. (H). 

The choice of the transition operator O in Sec. [VB 



implies that the second moments of the blocks in (|25 ) are 
(for large Nj) 



0^j)=N^/N' 



9h i 



0,1 



and 



NqNi 

iV2 



5051- 



(60) 



(61) 



Inserting this and the definition (|4^) of the functions 
Kj{ri) into ([5^) leads to 



(56) -r^ 



00 



[(1 - g.rff + r^^gf] (3] + 2if{l - gir^^)g^l3l 



= [(1 - g^iff + 77^50] /3| + 2??'(i - goif)goPl 

Wi = {^- 2gogiv' WPi + (1 - rf + 2gogiv')(3'c ■ 

(62) 

Let us specialize the se res ults to the case of the numer- 
ical simulation in Sec. IV B and see whether both are in 
qualitative agreement. 

The numerical simulation has been performed with 
iVo = iVi, whence, 



90 



1 

.91 = 2- 



(63) 



Since the total strength of the operator O is immaterial 
— see the discussion of Eq. ( p9[ ) — one can set /3j — 1 
which makes Pc = P- AH this turns Eqs. (^2|) into the 
more transparent expressions 



C2 - 1 - 
"-^00 ^ ^ 



1 z 
2^ 



(l-/?2) 



2 _ ri2 



'-^01 



*-^10 



1 - 

^01 



2^ 



(1-/3^) 



(64) 



If /? = 1, the operator (|2^) completely mixes config- 
urations with different isospin. In this case, Eqs. ( |64| ) 
yield 



'^00 



*-^ii — *-^oi ^ *-^io ^ ^ • 
The distribution ( p9|) of Wi / then becomes 

P{W.j)^G{W.,f,l), 



(65) 



(66) 



i.e. a Gaussian for all rj and, hence, for all a. This qualita- 
tively agrees with the results of the numerical simulation, 
see Tab. |[ 

If a = 1 (and Nq — Ni), the Hamiltonian ( |2^ ) is a full 
GOE matrix and its eigenstates have no isospin symme- 
try. According to Eq. (pO|), one then has r]^ = 1, Eqs. ( |6^ ) 
yield 



*-^oo — *-^ii — *-^oi ^ *-^io 



(67) 
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and the distribution (B9h becomes 



P(W,/) = G 



l/2\ 



(68) 



i.e. a Gaussian for all (3. Again this qualitatively agrees 
with the results of the numerical simulation, see Tab. H 



If neither the Hamiltonian (g3|) nor the operator (|2S 
mixes isospin, i.e. if (3 and 77 are zero, Eqs. (p4h give 



01, = 01, = 0. 



(69) 



In this limit, the distribution (M) takes the form 



W/) = i^{G{W,i,l) + W/))- (70) 



This says that half of the transition matrix elements is 
zero — namely the ones, where i and / have different 
isospin. The other half is Gaussian distributed. This 
result is expected and agrees with the numerical simula- 
tion, see Tab. || and Fig. ^A. 

In Fig. H the distribution of Wif is displayed for 
the mixing parameters a = 0.003 and (3 = 0.06. The 
hist ogram is the result of the numerical simulation in 
Sec. |IV I3| . The curve is the qualitative model (|5^). The 
parameters follow from Eq. ( ^ ) which yields = 0.0189 
and from Eqs. (p3) which give 



O?! = 0.981 , 
= Of^ = 0.0223 . 



(71) 



This illustrates that the qualitative model reproduces the 
essential features of the numerical data although it does 
not describe them quantitatively. 
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FIG. 6. A distribution of the transition matrix elements 
Wif. The histogram is due to the numerical simulation de- 
scribed in Sec. IVB. The curve is the qualitative model (59) 
with the parameters (71). 



VI. INTERPRETATION OF THE 
EXPERIMENTAL DATA 



In Sec. VIA, the experimental data described in Sec. || 
are compared to the predicti on of the qualitative solu- 
tion (H) of the model of Sec. |IV4 In Sec. |VI B| , we dis- 
cuss the effect that symmetries other than isospin have 
on the data. Parity conservation is used as an example. 



A. Isospin breaking in ^''Al 



The results of Sees. IV and suggest that the distri- 
bution of the reduced electromagnetic transition proba- 
bilities — displayed in Fig. |l| — does carry a signature 
of isospin breaking. In order to verify this, we compare 
in the present section the random matrix model to the 
experimental data. 

The parameters of the random matrix model are cho- 
sen as follows. The dimensions A'o, Ni of the blocks in 
the Hamiltonian (E3) are not equal. Rather one has 
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No/Ni 



(72) 



FIG. 7. The relation between a and 77^ for iVo / iVi. 



which approximately corresponds to the ratio of the den- 



sities of states with T — and T — 1 
entails 



90 = ^, 



1 



'Al. This 



(73) 



The parameter which quantifies isospin breaking in the 
Hamiltonian was set to a = 0.028. This value has been 
taken from the analysis of the eigenvalue statistics in 
Ref. 1^. We must convert it to the parameter 77 of the 
qualitative model of Sec. However, it is so far not clear 
whether Eq. ( ^3| ) applies because that function is based 
on a numerical simulation with A'o — Ni which is differ- 
ent from ([7^). We have therefore performed another nu- 
merical simulation of the distribution of the eigenvector 
coefficients a. This simulation was done with A'o = 150 
and A^i = 50. For each a, the function (^) was again 
compared to the numerical distribution of a and the op- 
timum 77 was determined. This yielded the crosses on 
Fig. 0. It turns out that the function (^ ) again repro- 
duces the numerical data if one considers only the range 
oi X = a/D ^ 1. Note that in the case at hand, the level 
distance is 

D = 0.192 . (74) 

For X — a/ D > 1 the function ( |4^ ) does not reproduce 
the data of Fig. Q because the Hamiltonian ( p3|) does not 
become a GOE matrix for any value of a ii Nq Ni. Due 
to Eq. (^4|) the numerical constant in ( |4^ ) is no longer 
equal to the level distance. 
One finds 

T]^ = 0.163 (75) 

for the present case of ^^Al. 

Due to the fact that ^^Al is a mirror nucleus, i.e. it has 
the isospin projection Tz = 0, the transition operator is 
uniquely determined by its isospin structure. 



If O is isoscalar it can by definition not connect states 
with different isospin. The parameter f3c in ( |2^ ) is then 
zero and, of course, /3 = 0. If the operator is isovecto- 
rial it cannot connect states with the same isospin. This 
is obvious for states having T = 0. For states having 
T = 1, this follows from the symmetry properties of O in 
isospin space. See Chap. 9 of Ref. specially the dis- 
cussion following Eq. (9.107). The parameter /3/ in ( |2^ ) 
is then zero. We evaluate the distribution (^) for these 
two cases. 

Let the operator O be isoscalar, whence f3c = 0. We 
set /?/ = 1, make use of Eqs. ( |73| , [75| ), and obtain from 
Eqs. (H) 



O^o = 0.922 
Of^ = 0.785 



Ol„ 0.153 



(76) 



This gives the distribution (|5S|) of the isoscalar matrix 
elements Wif. It is transformed to the variable 



(77) 



The result is the dot-dashed curve on Fig. g. 

Now let the operator O be isovectorial, whence /?/ = 0. 
We choose (3c = ^, use Eqs. (75j7^) and obtain from 
Eq. (p^) 



Ogo = 0.0782 



0.235 
0.847 



(78) 
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This gives the distribution ( p9[ ) of the isovectorial ma- 
trix elements Wjf. Again this distribution is transformed 
to the variable ([77[). The result is the dashed curve on 
Fig. |. 

One sees that both, the distribution with t — IV and 
the one with r = IS, are shifted with respect to the 
Porter-Thomas distribution (|^). This is a consequence 
of the appearance of Gaussian distributions with differ- 
ent widths in Eq. (|59|). The shoulder on the r.h.s. of 
the isovector distribution in Fig. |8| even hints to the ap- 
pearance of two peaks as in the numerical simulations on 
Fig. p3 and C. However, in ^^Al, isospin mixing is too 
large (namely A = 0.146) as to let the peaks separate. 
Remember that it is the class of small matrix elements 
that produces the peak which is centered at negative val- 
ues of z. The class of large matrix elements leads to 
a peak which is centered close to z = 0. The smaller 
the relative strength of the small matrix elements is, the 
more clearly the peaks separate. 
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FIG. 8. The distribution of reduced transition probabili- 
ties in ^®A1. The experimental histogram is the same as in 
Fig. 1. The full curve is the Porter- Thomas distribution (9). 
The remaining curves have been obtained with the qualita- 
tive model (59). The dashed curve — labeled r — IV — 
has been obtained with the isovector transition operator; the 
dot-dashed curve — labeled r — IS — is from the isoscalar 
transition operator; the dotted curve — labeled — is the 
weighted sum of both. 



There is a factor other than a (or ry^) which determines 
the shape of the distributions on Fig. |8[ For the isoscalar 
problem the class of the large matrix elements is the one 
with variance Oqq, see Eqs. (vW. It has the largest weight 
in the distribution, namely g^see (^,^^. Therefore the 
distribution is shifted only slightly to negative values of 
z. For the isovector problem the class of the small ma- 




trix elements is the one with variance Oqq, see Eqs. (78), 



which has the largest weight. Therefore the distribution 
is shifted more strongly to negative values of z. This 
makes the isovector matrix elements more sensitive to 
isospin breaking than the isoscalar matrix elements. 

The distributions for the isovector and isoscalar cases 
have been summed with the weights 530 and 343, re- 
spectively. These are the numbers of observed transi- 
tions, see Ref . Q . The result is the dotted curve on 
Fig. It is shifted with respect to the Porter-Thomas 
distribution by about the same amount as the experimen- 
tal data. However, the data are more strongly peaked 
around z « —0.5 than the present theory. This discrep- 
ancy must be attributed — at least partially — to the 



way the detection thresholds have been taken into ac- 
count, see Eq. (|l0|): The percentage of undetected tran- 
sitions may have been underestimated because its esti- 
mation was based on the Porter-Thomas distribution. 
The distributions that include isospin violation predict 
a larger probability of events outside the range of the 
experimental data than does the Porter-Thomas distri- 
bution. 

We conclude that the present investigation supports 
the conjecture of Ref. j|] which attributes the discrep- 
ancy between the Porter-Thomas distribution and the 
experiment to partial isospin breaking in ■^^Al. 



B. Other symmetries of the system: parity 

The nucleus ^^Al has more than only isospin symmetry. 
Especially there are quantum numbers which are exactly 
conserved or nearly so. For definiteness let us discuss par- 
ity conservation. Should it be included into all the above 
considerations that have focussed on isospin conservation 
alone? 



The matrix model of Sec. [V A can be extended to take 



parity into account. If the system has two possible values 
of isospin T = 0, 1 (isospin being broken to some extent) 
and parity tt — ±1 (parity being a very well conserved 
symmetry), the Hamiltonian H as well as the operator 
O can be written in the following block form 



(79) 



The four larger blocks are related to parity symmetry: 
The two blocks on the diagonal correspond to transitions 
between positive and negative parity configurations, re- 
spectively. The off-diagonal blocks correspond to parity 
violating transitions. They are zero in the present con- 
text. 

Each of the larger blocks is broken down into four 
smaller blocks related to isospin symmetry. They cor- 
respond — in a obvious way — to isospin conserving and 
isospin violating transitions. 

The "color" of the blocks schematically shows the 
strength of the matrix elements. The black blocks con- 
tain large matrix elements — allowed by both isospin and 
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parity conservation. The dotted blocks contain small ma- 
trix elements — allowed by parity conservation but for- 
bidden by isospin conservation. The white blocks contain 
zero matrix elements — forbidden by parity conservation. 

Since in the present context parity is considered an ex- 
act symmetry, the experimenters do not try to measure 
parity violating transitions. They assume them to be 
zero and — this is the important point — they do not 
even list them on their records. Hence, the experimental 
data set consists of the matrix elements in the two larger 
diagonal blocks. They are thus equivalent to the data 
that we have assumed in all foregoing considerations. 

Hence, the existence of exactly conserved quantities 
in addition to the partially conserved isospin quantum 
number is irrelevant if they are not measured. The exis- 
tence of partially conserved quantum numbers other than 
isospin would, however, alter the above results. 



VII. SUMMARY AND CONCLUSION 

In the present study, the effect of a partially bro- 
ken symmetry has been investigated for the distributions 
of wave function components and transition probabili- 
ties. This has been done by means of a random matrix 
model similar to the one introduced by Rosenzweig and 
Porter ||^. However, in the present model not only the 
Hamiltonian but also the transition operator may break 
the symmetry. This was dictated by the experimental 
material at hand and the symmetry that we had in mind. 
The data consist of a set of reduced electromagnetic tran- 
sition probabilities in the nucleus ^^Al. The symmetry 
is isospin conservation. Eletromagnetic transition opera- 
tors are in general not isospin conserving. 

A numerical simulation has shown that neither the 
components of the eigenfunctions nor the transition ma- 
trix elements follow Gaussian distributions. This is coun- 
terintuitive since the central limit theorem seems to re- 
quire Gaussian distributions if the wave functions are 
"sufficiently complicated" . Instead, the numerical sim- 
ulation suggests that the ensembles of both, the compo- 
nents of the eigenfunctions and the transition matrix el- 
ements, comprise at least one class of small elements and 
one class of large elements. The physical origin of these 
classes is clear: The small elements are isospin violating; 
the large elements are isospin allowed. At best, each of 
the classes has a Gaussian distribution. The probability 
density of the whole ensemble is then a superposition of 
Gaussians with different parameters. 

This idea has been worked out in detail and results in 
a qualitative analytical evaluation of the random matrix 
model. It yields simple formulae for the parameters of 
the Gaussians. Its practical use is to circumvent easily 
the demanding numerical simulations. Conceptually, it 
offers to see in detail, why the premises of the central 
limit theorem are not given: The elements in the dif- 
ferent classes depend in different ways on the symmetry 



breaking parameter a as well as the dimension N of the 
space. It is well known that in the limit of TV — > oo an 
ever smaller value of a suffices to make the above classes 
indistinguishable and thus to produce simple Gaussian 
distributions of the eigenvector components and the tran- 
sition matrix elements. However, this is what happens if 
N is taken to infinity such that the parameter A = a/D 
goes to infinity. If N is finite or if it is taken to infinity 
together with a ^ such that A remains finite, then the 
above mentioned classes conserve their identity. It is this 
situation which the present paper deals with. 

The present results thus support a conjecture formu- 
lated in Ref. Q on the basis of the data collected from 
^^Al: The authors conjectured that the distribution of 
the electromagnetic transition probabilities carries the 
signature of partial isospin breaking. The present theory 
offers a qualitative understanding of the data. Whether 
it can quantitatively reproduce the data remains to be 
seen. We have not attempted to fit the theory to the 
data. 

A caveat seems necessary: The random matrix 
model ( p3[ ) applies to any partially broken quantum num- 
ber — it is not specific for isospin. Certainly, the avail- 
able information on isospin breaking in ^^Al requires that 
the distribution of the transition probabilities differs from 
a Porter-Thomas distribution to the extent given by the 
dotted curve on Fig. ^ — but the data may in princi- 
ple and in addition carry the signature of other partially 
conserved quantum numbers. These may not be as fun- 
damental as isospin but rather given by the specific form 
of the effective nucleon-nuclcon interaction in the sd shell. 

Very recently, an analytical calculation, based on the 
supersymmetry method, of the wave function statistics in 
coupled chaotic systems was presented in Ref. . These 
results are valid for the case of broken time reversal in- 
variance and therefore, unfortunately, are not applicable 
to the case of conserved time reversal invariance which is 
discussed in the present paper. 
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